00۵ 0 ۷۲۷۵۵۲ 2۵0 0 8 


کع صصن۵ا۰ ۵۹۰/۸19 


6 62020[ 
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0 )606۵و و۲۱۵۲ صرز م۱۷۵0 ۲۷۷۸۵ ۵۲ ماع(مصه هام۲۵ 20 از تاتعصهو 
۷۷۵۱60 عمصنهامهم]۷ ۵ ۲۵ ممت)ما بصن یمام عمط 
0۵660606 ۲۲۵۲ زو اه 0۶ تا ععهن)) 


۱ 


:6 318 ماه 0 ۲۱۵۱۷ 16-10-01 :16061۷60 
مصح ۵25101۷10۲ .2022 .1 نطاج۳ مه م۷۲ تصوزفک ریق 24ظ2مصیتم۲ 04-01-02 :1۹6۷1560 
کمطحومطظ مه )صمحصتلعو رما ما ۵0۵1 5۷/۸۲ ۶ه دوب امه اصتدانهعملا 10-04-02 :۸۵۰۵۵0۲۵0 
۷۵ زه۵۲ک۱ ۵۶ تاو عققن)) 04عطوع20 ۷۷ ویامصتماصتم۷ و 10۲ صمت)ها نار 21-07-02 :عصنلمن فاطقاند۸۱ 
۳۵۳۹1۵۴ 12) ,167-183 :(36)02 الم مضه عهاج/۷۷ ۵۶ امصننمد. .,فصمصصطء‌امی 
0اعحتاوطاه طفناعطظ 
00 :]10 


۱۱۹۵0۵0 

۲ هط ۵۶ عصم و1 و48مطوععمنه فبامصتهاصنمص صا امصند ممهسه ۵ وتلمتان ‏ مه واتاصمتن ۵ ممتاه‌ااصل 
۵۵ ۵۵۱۵ مه طام‌مته۵م۵) آمیاوباصنا قق طمناه رقععنا2ع؟ مبامتصت وا ما بل عصتاممصظط ما وملمم] مصتصما مه 
2۵۲ ۶ ولورلقمه اصته)عمصن مطا و فاصعصصطمق میاه عصتاممصه ۵ فامعمو2 16996-1000۵ مطا ۵۶ مق .2۲0689965 
امن بصماتمموصا مط ۵۶ موم ود فنامطامومطاظ تمامصصهتهم فتمطام‌ومهام هار۷ عمط یهام موه رفهمتامت0ع۳ج اتلمتان 
40ج عصعع‌تاه ما فمممممممطی ممتهمتمنمتانه مفاهه صوع. وعم‌تیاموع1 تم۳۷/2 18 12076286 فلا 24 رتعا2/ طا ۱۷۵۲1۵۵0165 
عصناه0مصظ من تمتجه قطا مج عطا صرح ععاعصصدتهم فنمطمومطام فطا ۵۶ عمعصه)زمم‌صا فطع م6 ما1۲ .عصصقل ۵۶ ۲۵96۲۷01۲5 
5 ۱۵06160 معط فقظ فباعمطمومطام م0 رط‌تهعوع: ونط صا رد 0و 2ععه اصمحصطماق زدتقکا عطا مر ۵0ع«مامصصه صععه )مظ فقط 
۲ فطا تاصع10 ما فصصتع لاه فتط1. ,ععتهطهون اصمصصتعی هه م1 مطا ط عجصماه متمصصهتهم اتمتان ماه 2 
40ج ماهنل۷2 رماه:هه 0صه معتمطمون فیمجامومصام ناما فصه رتصمحصتلعه سمل مه مق فص ۵۶ وتعامصصهتهه 9699101۷70 
کنامصتق)صیامحط 2 طاً وماطاهتنه۷ معط عفقط مصتامن0عنم صا (ع0مصص 5۱۷/۸۲ فطا ۵۶ باصتدانهعصنا متعاعصصه‌تهم فص 22728 
۶ 086 و1 من بع0۵ صدوملاظ ۵۶ دون تع۴ زجتهک! معط ۵۶ م2 اممصصهه‌هه معط مه لاه عفقع م1 امه 
۷ ۱۳۲۷۵۲ نفک مطا عصماه فمتاماه ماه لهعتاتن ۳0 تچ م۲۵2۲ بصقا صا فلمطوماهه فنامصتهصتمقظر فط 
۷ ۱۵ 009۲۱0۱۵0 و1 کقطا اععزم0 096صرتناج- تج 2 وا رژهتفک) صفل اتمه لصمه صحمااظ قمع صقل هتشر 
صحوعاظ م1 .رتمک ۵۶ وطاتناطانه مطا صا متبقلنم‌تهه 4ص ممتاهع1۳1 10۲ تماجه عقاوم 0ص صحعطه 1 مه ماه ومیل 
ونطا ۶ه عمصها۲محرحصصا فطع من بصحتطاه] صا ماه مصتلصتیل مه تمتممه مد «آم‌مرناه ۵۶ نصتمم تمتاصمووع ما مفله و1 ۵4 
0 24 ع۷تاهاتاصهتان ۲ع ع0مصه مامتطانلی 2 مصتلب۵عص روعوتا اصمعکنل ت10 ماه فصارام‌مناه صا ط210ع۲ 
11 فنطا صد دمتتقصمعو )ممجم‌مهصفصظ معتبظ راومه م) فصمتوتععل ۶۵۲ وتفوها مط ۵ صقع ۱/۵0۵۲ 0۴ ۱۷۵2۲1۵0165 


۵۵۵۹ 20 واهز۱۷]۸)۵۲ 

طونم معجه27 صه طاز من رصتفوظ محهلنظ که موم همه اممصصهم‌ناه ۷۵۲ زهتفک مطا مه تناو مهم مط 1 

۶ 2۳62 صح ۳1 صلجوه وتط1 .فصتقمتم]۱ ۵907 مط طاً 1062060 قهعته قتامصتعصتممه معط ۵۶ عمجم و رعتعاعطظ 2880 0۲ 
عط ۶ه طجامو مصح )عفع فطا ما .عمص1ه۳ج صهتهلصه۱۷]22 ۵۶ فاتدم 101065 بطاعمه فطا صا فعاعصمات۱! عتمباود 1076 
2۵۵۵ 132 ,0۲0۷۱۵۸۵۵ ۵۱00۲2 طا 10621060 و1 14 ۵۶ )دمص مه ممصز۵۲0۲ صهتطام 1 ۵۶ فاتوم وعلسامطا همه اممصطمهه 
7 0 2۷۵۲۵2۵ معط مضه ,باه۲۵۹06۵۲۷ رح 480 مضه ۰*6 12.1 ع27 صلفعه فنطا صا الملصنده 4ص متنه‌مم‌جها تمناصصه 
عط ما کصمعهز20 فصمتاهاه متاممصره 0 ماع اند صم-عمم۱ م1 بههعه ونطا صا ۷60موهاه و1 بقع مظ) عم دوم لمتمهاع 
,۸150 .۱۵061 قطا م 0۲000660ص1 قوب (رلقاما ما ویدع۷ 21) 2018 0۶ مه فطاٌ م6 1998 0۶ ععتصصتععه مطا حصم عته تاه 
۵ 2274 ,وتنام ممتهت۲2 تقامو روما حصیاحجقصظ مرج یتمه رالهگصنه؟ ووائ‌تحصط متام ۵ 0242 «اتدل 
۴ .۱0061 36 0 100000660 ۷۷۵۲۵ 2762 تاو فط ما فصمتاهای 2 60عناقجمصظ فعامصصه‌تهم 26160101021621 قه ٩0660‏ 


یمیمص ۵۶ صمصندم۵ 1 رتمووع۲۳۵ صهاوتومخه مصه مملهممتمصظ متممصمماطظ که ماه0۳20 ,۷۲,9۵ -2 4ص2 1 
۷۰ ,1۲21 یداه 1[ ,زالوتع نا وعتم۷۵0 )1ات2 1 رمصارممصاعمظ ام)صمصمم دص مرج زذ0 ۵۶ لمع روصتهعصنم۴۸ 
(۵۵02169.26.11) 201 رصح تاتقعهظ ممطه مصت‌ممصوعتزمن -۶) 

مد رصحتطه 1 راممصمم صظ ۵۶ اممحانه1۱۵۵ بمم0۶1 9011 4ص عماج ۷۷ عط) ۵۶ بمممصه] «انامعب[ -3 


۸ نشریه آب و خاک. جلد ۰۳۶ شماره ۲ خرداد - تیر ۱۴۰۱ 


۱ ۱0 
0 صصعل اماف‌اتصه هم 0معتمطونل ماه «مل «اتقه مممتهصمهملص ومطاج مطا م6 صمت0ع ظ هلق همم 
,10101۳02۵1108 کصمجدرمع 9212 0۳020 رحمتل20 1۱ ,مه فطع م6 1000006860 ۷۷۵۲۵ هل مطا ۵۶ فصمتاه»لمممصو همهم 
۵ ,0۲09۲19 لهصمنوع1 صا 0مویا وع۳01126ع1 متقطامومطن م6 ۳۵12660 مملافصمکصا 0ج فهمتهه ممت)اهعتت مصللتاامط1 
۶ مهن 69910۷۵ مطا مصتصماعل م۵ ۱2۵۵0 ۷۷2۵ 00طامصه فلورولمصه اد قمع ل92ماع م1 ,061مصظ مطا م4 )0۲۵968 
م6 2021776 مهصه مهتصاللهه 6۵ مهو ظنآن ۷۷۸۲ ص1 عون قه حصطات:0ع2۱ 9۱۲۳۲2 عطا رععمصصعطتن۳ ,1ع0مقظ فطه 
9590۵ 0۷ «اتهازهمصنا آنامدان فط فعتصهین صطاتمعل2 فنط1 .۵061 5۱۷۷۸۲ معط ۵۶ بوامتماهعص متتامصصهته2 
۰ 0 اصتهازهعها 


تمه کز(۱ 20 و6اباععع1 

,(۳ ۳۱۵صله) ۲مامع همه س«مالهعوظ فمامصصهنهم مطا روت رافمه و تمصع ۵۶ عناوم عمط ما عصتل۲هععضظر 
6 ۳۱0۹ 66 ۷۵۲۵ (۳۱۸۳۲9) م22 ووقا ممتاهالممع۲ م2 رز( 0۷) ماک امن هه تلو و مصتممه]۱۷ 
۲ مهم (ولن) ناما مها اوع9 مظ1 .۵6و۲۵ رفتامطاندمدام هام 24 رصمصصتلعه سمل ما وتمامصطه22 
00 274 2)10هاناهه الب تمه مه فصمتاهاه اه مد 0عصتم‌هاه ممتافانصته ممدامومصای تماما مه رامه‌صتلعو رصن 
معط 211 ۵۶ حصحتع112 معط 0۶ ماصتمم لدعم معط ۵۶ ممتامانته .راه۲۵5۵6۵0۷ ,0.92 4صه ,0.56 ,0.76 6۵ ت60 ۷۷۵۲۵ 2۵۲1005 
م0 06 ]اه ,محتهعه ممتامت0ع7( [06مظ 0۵082880 هه «اصتههعصن 180762860 ۱/10 2980012160 0فل2 فه وعالاصهت۳ان 
۸ 1386 0۴ 60۶۵ 6271۷ مصه فمنله۷ )حمصصتلعی صرح ۶مصنار 0عتناعحعطه فطل ۵۶ 709۵ صقطا معمصط گم اصمصدرهمهام مطا 
۶ مبلد۷ امص معط حصطاترمعاه 5۱۳۳۲2 رها 0ع)همجهم ومتتهصنهه راصتهاتهمصرنا ممنام01عتم فط را ومیالهه فنمطام‌ومطم مها 
لته صا ومتتهجهعو فیط عصتافاله 10 ماقتتومتومه و9 صقه قومممم ممتاهعطانلهه مط هد 0عفیا فتمامصصح‌تهم معط 
اممصصتلعو لمح سمل ۵ فاصتمم علوعم مطا عصتاهابصدته ود عمج فطع ۵۶ وعمصلهه صتقمدر م1 .عاصمحصطمتق صتق/صتامظ 
6 ۱0 عبال 012117060وعع 0 مه عصتام0مصه ۵ طتااطاهنا مطا ,ععقفطهفنل تممصصتلعه قمع سمل ۵۶ موه مط 1 .عععقطهولل 
۷ ۱ 10616296 ۲0 فصتععه مامت رعصتنمو صا 1۷۵ ما مه «ماکصا تامصوعمصی مطا ۶ ممتامتلع؟0 ملحتتاهعه 01 م12 
صهه مصمتافاه لهملع۳۵۵۱۵0۲0۱0 طً 0202 عصتووتصه گم ممامصربنطر ماماهععلنموم هر ۱۷2۷۰ طا عصتمج علقعج فص وعطمهع1 200 
07 عصاوول رطمنعع۳ فنطا ما عاصتمم لقع مطا 2 ممتامتلع۳ج [06مصظر ماقتنهعع 0۲ عامق1 مطا. )۵166 باهتاهعرام 
عم 0۶ م0۵۵0 مط1 اصه‌تصونه مها فقظ موه لهملوما0۲0ممامص تفا ما جهن قاقل ناماما 
۵ 2 ۱۲9 0عه آل-و وه همه فاعهمصظ )منم هنن مصه اماصمصتهمره صقر نزن مق 5۷۷۸۲ 

۰ و2 ,0900272 )ممحصزلعو مطا 0۶ فاصتوص لقع مط) مصتامت0ع۲ج مد وومصاهع۷ مطا عصتامعله ص120 تمطا0هه 


«متعصلصم) 

فتمطامومطح لقاما 4ص ممصتلمو رما 0عموهان 0۶ دمص رقالتاوع فلورلقمه ماصتقاهعمه مطا ما عصتلتهععخظر 
ولا 16 ,صطاتیمع21 ٩۱۳۳۲۲2‏ عطع ها 0ماجعصمع ومتیه)مهه ممتاما0هتج واصتمه عم ما تاه م۵ ومیاله۲ ععتقطهولل 
عط نقطا قمع 16 ,۰۲۵۵۵۲01۲۵ ,10مطومنطه عصتاملمصه بمامدفتاهه مط افص فقط فماطهتنه۷ معط هب1۵ اممله‌تااههه 
6 صهه میاه تفص تتمطا مج دربتاه فتطا صا وقععمعم ممتاهتطاتله فمطا صر عون 0ص2 1060111160 فماممصههم 8699101۷70 
معط ۵۶ موم .عاممصصطماهه صتقاصمصظر مقلتصصته ‏ مضه معه تاه ونطا عم ومتنمجععو ماب عصتام0مصظ عم ع200۲001126 
مطا 101 «للملمعووه رفصمتاهاه ۱6160101022 طاً 024 عصتوعتمه ۵ بیط معقا م فهب بو عصعومتم مطا م۵ فصمتاهاتطتز 
6 ععصقطصصه بچحوه اممصه مه مه هه همزع010۱۵عامص 0متناقوعصه 0متنا0ع۳ مصللز ۳0ج رقتاطز. ماماهه۷ عتااهتهم‌طظما 

۰ 0621 26 656012117 رطمتاهلباطظله 


۷ 216۳ ۷۷ ,۱۲۵061 01500۵0 تم ر۳۲]2 ۱ رمتام(۲:۵0 ماهتا 16۲۱۷۵0۸۵۵۵۰ 


8 نشریه آب و خاک 
و ۷ 0( 


جلد ۲۶ شماره ۲ خرداد-تیر ۰۱۴۰۱ ص ۱۶۷-۱۸۲ 


تحلیل حساسیت و عدم قطعیت مدل ٩۷۷۸‏ در شبیه‌سازی دبی جریان. رسوب و فسفر برای 


یک حوضه آبریز کوهستانی (مطالعه موردی حوضه رودخانه کرج) 


۷ ۰ ۱۶۰ ۲ مه ۹ ۷۳ 
شاهین نوری نژاد - محمدمهدی رجبی "*- تورج فتحی 
تاریخ دریافت: ۱۴۰۰/۰۷/۲۴ 


تاریخ پذیرش: ۱۴۰۱/۰۱/۲۱ 


چکیده 


شبیه‌سازی کمیت و کیفیت رواناب‌های سطحی در مناطق کوهستانی» به دلیل ویژگی‌های منحصربه‌فرد اين‌گونه حوضه‌های آبریز مانند تأثیر 
قابل‌توجه برف بر رژیم جریان» تغییرات زیاد پارامترهای هواشناسی به‌واسطه تغییر ارتفاع, نرخ بالای فرسایش و کمبود اطلاعات ناشی از مشکلات 
دسترسی فیزیکی» یکی از موضوعات چالش‌برانگیز در حوزه مدل‌سازی حوضه‌های آبریز محسوب می‌شود. یکی از مواردی که کمتر در خصوص 
مدل‌سازی این گونه حوضه‌ها به آن پرداخته شده, تحلیل حساسیت و عدم قطعیت پیش‌بینی‌های مرتبط با کیفیت آب. به‌ویژه در ارتباط با پارامتر مهم 
فسفر است. هدف از انجام این تحقیق واسنجیء اعتبارسنجیء تحلیل حساسیت و عدم قطعیت پارامتری مدل ٩۱۷۸‏ در پیش‌بینی دبی جریان» رسوب 
و فسفر کل در یک حوضه آبریز کوهستانی است. مطالعه موردی در حوضه آبریز رودخانه کرج در بالادست آبگیر بیلقان انجام شده و الگوریتم 51/۳12 
برای انجام تحلیل عدم قطعیت مورد استفاده قرار گرفته است. بر اساس نتایج تحلیل حساسیت. پارامترهای ضریب جریان پایه آب زیرزمینی» ضریب 
زبری رواناب سطحی و آهنگ افزايش بارش به ترتیب حساس‌ترین پارامترهای مدل نسبت به دبی جریان» رسوب و فسفر کل بودند. بهترین ضرایب 
نش-ساتکلیف (5) مربوط به شبیه‌سازی رواناب» رسوب و فسفر کل در همه ایستگاه‌ها و در مجموع دوره‌های واسنجی و اعتبارسنجی به ترتیب برابر 
با ۰۰/۷۶ ۰/۵۶ و ۰/۹۲ بدست آمد. همچنین شبیه‌سازی نقاط اوج نمودار هر سه کمیت مذکور با افزايش عدم قطعیت و کاهش دقت پبش‌بینی مدل 
همراه بوده» اما با توجه به قرارگیری بالای ۷۰ درصد مقادیر اندازه‌گیری شده رواناب و رسوب و نزدیک به ۶۰ درصد مقادیر اندازه‌گیری شده فسفر کل در 
بازه عدم قطعیت پیش‌بینی تولید شده توسط الگوریتم 51/۳12 محدوده پارامترهای مورد استفاده در واسنجی مدل می‌تواند الگوی مناسبی برای 


شبیه‌سازی سناریوهای آتی در حوضه‌های آبریز کوهستانی مشابه باشد. 


واژه‌های کلیدی: پیش‌بینی. کیفیت آب. مدل نیمه‌توزیعی» واسنجی, 51/۳12 


مقدمه 


در دو دهه اخیر مدل‌های نیمه‌توزیعی هیدرولوژی و کیفیت آب به 


تربیت مدرس» تهران» ایران 
(#:- نویسنده مسئول: ۴ 6۵ ماه زهتمندط. :اتف 
محیط زیست» تهران» ایران 

10]: 0 


طور وسیع برای تصمیم‌سازی جهت مدیریت حوضه‌های آبریز و 
پیش‌بینی آثر وقایع طبیعی و فالیت‌های انسانی بر کمیت و کیفیست 
آب به کار گرفته شده‌اند. یکی از پرکاربردترین مدل‌ها در این زمینه» 
مدل 5۱۷۸۲ است (1998 .۵1 » ۸01۵). نتایج تحقیقات پیشین 
حاکی از عملکرد مطلوب مدل 5197۸7 در شبیه‌سازی فرایندهای 
مختلف از جمله انتقال رواناب سطحیء رسوب و پارامترهای کیفی در 
مقیاس‌های مختلف بوده است (2007 .۵1 61 0۲و۸0 : 
6 .0 6 12۳002 :2013 .]6 6 ۱ در مدل 5۷۷۸۲ 
دا گناد با رامق ورخدی موه درد که سای از آنخ‌ها با اساسا جه 
شکل مستقیم قابل اندازه‌گیری نیستند» و يا اطلاعات مناسب برای 
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توصیف کامل تفییر آن‌ها در سطح حوضه آبریز عموماً وجود ندارد 
(2010 ,.]۵ ۶۶ ط۲216012206) و لذا آنچه که به عنوان ورودی مدل 
به کار گرفته می‌شود. تنها تخمینی از مقادیر این گونه پارامترهاست؛ 
بنابراین تحلیل عدم قطعیت مدل 917۸ به‌منظور ارزیایی کارایی 
نتایج آن در شبیه‌سازی‌ها امری اجتناب‌ناپذیر است. 

شبیه‌سازی حوضه‌های آبریز کوهستانی با استفاده از مدل‌های 
نیمه‌توزیعی مانند ۹۷۷۸۲ به چند دلیل یکی از موضوعات 
چالش‌برانگیز در حوزه مدل‌سازی محسوب می‌شود ( ,./4 6۶ ۷۱۷1۲011 
7 228 ۳۱2829 :2015 .]6۵ 6۶ ۷۷۵ :2014 .۵ 61 ۱۵0۲ : 
7 ,ع۳۳202020). نخست اینکه تغییرات مکانی پارامترهای 
هواشناسی مانند بارش و دماء به دلیل تغییر زیاد ارتفاع در سطح 
حوضه‌های آبریز کوهستانی» اغلب بسیار بیشتر از دشت‌ها و مناطق 
کم‌ارتفاع است که اين موضوع تأثیرات قابل توجه در پیچیده‌تر شدن 
فرایندهای هیدرولوژیکی دارد (2010 ,اه 6 نصقاان ملق : 
6 :00۵1 0ص2 طعصت۹ :2011 ,مصطه 0صه). دوم اینکه ذوب 
برف در ارتفاعات منشاً مهم ایجاد رواناب‌های فصلی است و تأثیر 
زیادی در رژیم فصلی جریان آب دارد (2008 ,عحصهئل ۷۷ 4ص ن : 
5 , .4 61 01۷5505)). این موضوع نیز شبیه‌سازی با استفاده از 
مدل‌های رایج هیدرولوژیکی را با محدودیت مواجه می‌سازد 
(2013 .۵7 6۶ 2020). سوم اینکه مناطق کوهستانی نسبت به 
سایر نواحی با خطر فرسایش جدی‌تری مواجه هستند 
(2012 ,اه 6 تحصصنامصصهوه:۳) و سرعت فرسایش و انتقال 
رسوب در رودخانه‌های کوهستانی بسیار بالاست که زمینه‌ساز افزایش 
غلظت آلاینده‌هایی مانند فسفر می‌باشد (2011 .۸1 2 2۷۲۵/1 : 
3 , .0 2 ۷۷۷). کمبود اطلاعات ناشی از دشواری دسترسی 
فیزیکی نیز چهارمین عامل چالش‌برانگیز در مدل‌سازی حوضه‌های 
آبریز کوهستانی بشمار می‌رود. 

به واسطه اين دلایل» پژوهش‌های متعددی در سال‌های اخیر به 
طور خاص به بحث مدل‌سازی نیمه‌توزیعی حوضه‌های آبریز 
کوهستانی پرداخته‌اند (2018 .۵ 6 ع2 ۷۷ ,.۵1 61 تلحصوتب0ع۸ : 
1 , .۵ ۵ عط2ط2 :2020 .۵ 1 16088 :2019) و به کمک 
این پژوهش‌هاء از یک سو دانش مدل‌سازی حوضه‌های آبریز 
کوهستانی پیشرفت نموده و از سوی دیگر محدودیت‌های مدل‌های 
رایج مانند ٩۱۷۸۲‏ در شبیه‌سازی این‌گونه حوضه‌های آبریز نمایان 
شده است. یکی از خلاء‌هایی که در مطالعات قبلی وجود دارده عدم 
توجه به بحث تحلیل حساسیت و تحلیل عدم قطعیت پارامتر کیفی 
"فسفر" است. این پارامتر از آن جهت که اغلب عامل محدودکننده 
تغذیه گرائی در منابع آب شیرین بشمار می‌رود. دارای اهمیت ویژه در 
مباحث کیفیت آب است (2008 ,7ع2:0601))» و در حوضه‌های آبریز 
کوهستانی با نرخ فرسایش باله به میزان زیاد وارد مابع آب سطحی 
می‌شود. حوضه آبریز رودخانه کرج یک نمونه از این دست حوضه‌های 


آبریز کوهستانی است. میزان فسفر در آب رودخانه کرج به دلیل 
پوشش گیاهی نسبتا فقیر و نرخ فرسایش بالای حوضه آبریز 
به خصوص در فصل بهار قابل توجه است و این موضوع سبب شده تا 
آب سد امیرکبیر در پایین‌دست این حوضه آبریز» در مصرض 
تغذیهگرائی قرار گیرد (1998 ,.۸1 2 ۸۳014 ,.اه ۶ تتتعط 

:2018 از طرف دیگره حوضه آیریز رودخانه کرج نقش بسیار مهمی 
در تأمین آب شرب پایتخت ایران (تهران) ایفا می‌کند و لذا مدیریت 
بار آلودگی فسفر در این حوضه از اهمیت فوق‌العاده برخوردار است. 

در اين راستاء مطالعه حاضر به شبیه‌سازی دبی» رسوب و فسفر 

کل در حوضه آبریز رودخانه کرج با استفاده از مدل نیمه‌توزیعی 
۲ می‌پردازد. تاکنون مطالعات متعددی پیرامون ارزیابی کارایی 
مدل 5۱۷۸ در شبیه‌سازی‌های هیدرولوژیکی و مدل‌سازی رسوب 
در حوضه‌های آبریز کوهستانی انجام شده است. برای مثال فلین و 
ون لیو (2011 ,167 ۷2۵ 4هد صصرا۳) از ٩۲۷۸۲‏ برای 
شبیه‌سازی دبی جریان و رسوب در حوضه کوهستانی برف‌خیز لامار! 
در آمریکا استفاده کردند. یافته‌های آن‌ها نشان می‌دهد که 5۷۷۸۲ 
برای شبیه‌سازی هر دو متفیر در مناطق کوهستانی و تحت سلطه 
ذوب برف مناسب است. آندریاناکی و همکاران ( .1 4 تلتصهتته۸ 
9 اذعان داشتند که مدل 5۷۷۸۲ در منطقه مورد مطالعه آن‌ها 
زمان‌بندی ذوب برف را به طور دقیق پیش‌بینی کرده اما نتوانسته اوج 
جریان را برای سال‌های خاص با دقت مناسب شبیه‌سازی کند. از 
معدود مطالعاتی که به شبیه‌سازی فسفر در مناطق کوهستانی 
پرداخته‌اند می‌توان به مطالعه لیو و همکاران (2016 ,.۵1 6 11) در 
حوضه آبریز رودخانه شیانگشی " اشاره کرد که توانستند با استفاده از 
7 پیش‌بینی‌های خوبی از فسفر ارائه دهند اما در نتایج آن‌ها 
تحلیل عدم قطعیت فسفر منعکس نشده است. در این تحقیق برای 
اولین بار در یک حوضه آبریز کاملاً کوهستانی که بیشتر سطح آن را 
مرتع با پوشش گیاهی کم تا متوسط در بر گرفته و فعالیت زراعی در 
آن به نسبت مناطق پست بسیار محدود است. تحلیل حساسیت 
جهانی " و تحلیل عدم قطعیت پارامتری مدل 5۷۷۸۲ در شبیه‌سازی 
ماهانه فسفر کل انجام شده و تأثیرگذارترین پارامترهای مدل در 
شبیه‌سازی فسفر کل (به همراه دبی و رسوب) معرفی شده‌اند. لازم به 
ذکر است که مدل‌سازی دبی و رسوب پیش‌نیاز مدل‌سازی فسفر کل 
تب بافتاه 
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مواد و روش‌ها 
منطقه مورد مطالعه 

منطقه مورد مطالعه در اين تحقیق یعنی حوضه آبریز رودخانه 
کرج تا محل آبگیر بیلقان» با ارتفاعات ۱۳۶۴ تا ۴۳۶۹ متر و متوسط 
ارتفاع ۲۸۸۰ متر جزء مناطق کامللاً کوهستانی است که در سلسله 
جبال البرز واقع شده است. این حوضه با مساحت ۱۰۷۶ کیلومتر مربع 
در شمال خود بخش‌هایی از استان مازندران» در شرق و جنوب حوضه 
بخش‌هایی از استان تهران را در بر گرفته و بخش اعظم آن در استان 
البرز قرار دارد. حدود ۶۰ درصد از این منطقه دارای ارتفاعات ۲۰۰۰ تا 
۰ متر از سطح دریا می‌باشد. میانگین دمای سالانه و بارندگی در 
این حوضه به ترتیب ۱۲/۱ درجه سانتی‌گراد و ۴۸۰ میلی‌متر است و 
میانگین ۱۱۷ روز یخ‌بندان در طول سال در این منطقه مشاهده 
می‌شود. در مسیر رودخانه کرج دو سازه آبی بسیار مهم یعنی سد 
امیرکییر و آبگیر بیلقان قرار دارد. سد امیرکبیر (کرج) یک طرح 
چندمنظوره می‌باشد که از مهم‌ترین اهداف آن تأمین آب شرب تهران 
و تنظیم آب برای مصارف آبیاری و کشاورزی اراضی حومه کرج 
می‌باشد. همچنین آبگیر بیلقان مهم‌ترین نقطه تأمین و انتقال آب 
شرب تهران به حساب می‌آید. در شکل ۱ موقعیت جغرافیایی منطقه 
مورد مطالعه مشخص شده است. اراضی واقع در محدوده مورد مطالعه 
شامل کاربری‌های پیکره‌های آبی اراضی باغی " و مرتع است که از 
نواع مرتع سه دسته مرتع با پوشش گیاهی ضعیف " متوسط " و خوب؟ 
در این حوضه وجود دارد. بیشترین کاربری اراضی در محدوده مربوط 
به کاربری مرتع با پوشش گیاهی متوسط است که مساحتی در حدود 
۷ کیلومتر مربع و معادل ۸۶ درصد سطح کل محدوده مطالعاتی را 
در بر می‌گیرد (شکل ۲). بر اساس نقشه تهیه شده از نقشه جهانی 
خاک فائو » منطقه مطالعاتی دارای سه نوع خاک می‌باشد که در شکل 
۳ قابل مشاهده است. 


مدل 5۲۷۷۸۲ 


۲ یک مدل هیدرولوژیکی» پیوسته زمانی» نیمه‌توزیعی و با 
پایه فیزیکی است که توسط سازمان تحقیقات کشاورزی ایالات 
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متحده آمریکا توسعه داده شده است (1998 .۸ ۶ ۵010 این 
مدل مفهومی به‌منظور پیش‌بینی تآثیر روش‌های مختلف مدیریتی بر 
مقادیر آب» رسوب و مواد شیمیایی کشاورزی در سطح حوضه‌های 
بزرگ و پیچیده با خاک» کاربری اراضی و شرایط مختلف مدیریتی و 
مورفولوژیکی در درازمدت طراحی شده است. ٩۷۷۸۲‏ قابل اتصال به 
نرم‌افزارهای 15) بوده. لذا حجم وسیعی از داده‌ها به‌صورت لایه‌های 
اطلاعاتی مانند نقشه کاربری اراضیء نقشه خاک و نقشه ارتفاعات 
رقومی " (/0۳) با استفاده از نرم‌افزار واسط 015 توسط مدل بکار 
می‌رود. مدل 5۷۷۸۲ برای شبیه‌سازی حوضه ابتدا آن را به چند 
زیرحوضه تقسیم می‌کند. هر زیرحوضه نیز به چند قسمت دیگر به نام 
واحد پاسخ هیدرولوژیکی" ([11181) تقسیم می‌شود که این واحدها 
قسمت های درون هر زیرحوضه هستند که دارای نوع خاک» کاربری 
اراضی و شیب یکسانی هستنن. لین ترئییب از نظر مکتانی: 111007 
کوچک‌ترین واحد محاسباتی مدل 517۸7 است. برای شبیه‌سازی 
فراین‌دهای هی _درلوژیکی و کیفیت آب 5۱۷۸1 از معادلات و 
روش‌های متعددی استفاده می‌کند که به‌تفصیل در منابع متعدد به 
آن‌ها اشاره شده است (1975 رعصعا۷۷11 :1977 ,0امحعدظ : 
72 .۸ 6 ۸۳9010 :2011 .]۵ ۶۶ ط0فالع : 2005 ۱۱۵ 


راه‌اندازی و اجرای مدل 

در شروع کار با مدلء ابتدا نقشه 9 که قبلاً تهیه شده بود» به 
مدل معرفی شد. پس از ترسیم شبکه آبراهه‌های محدوده مطالعاتی بر 
اساس 17۷( معرفی شده حوضه آبریز توسط مدل ترسیم شده و 
زیرحوضه‌های تشکیل‌دهنده آن نیز مشخص گردید. با درنظرگرفتن 
حداقل سطح زهکش ۷۰۰ هکتار برای ترسیم زیرحوضه‌ها» تعداد ٩۱‏ 
زیرحوضه ترسیم شد. در مرحله بعد با استفاده از نقشه‌های رستری 
مربوط به خاک کاربری اراضی و تعیین طبقات شیب واحدهای پاسخ 
هیدرولوژیکی تعریف گردید. برای این کار دو نقشه اول از قبل آماده 
شده بودند و طبقات شیب نیز در همین مرحله با استفاده از مدل 
و در محبط ۸۰015 تهیه شدند. با توجه به تنوع بالای 
شیب در منطقه, پنج طبقه شیب ۰-۱۰ ۱۰-۲۰ ۲۰-۲۵ ۳۵-۵۰ و 
بالاتر از ۵۰ درصد در مدل تعریف گردید. همچنین با توجه به آهمیست 
تعریف بازه‌های ارتفاعی به‌منظور بهبود شبیه‌سازی در نواحی 
کوهستانی (2013 ,.۵1 6 صححصطف؟ :2020 ,.۵ ۵1 202ظ ) تعداد 
پنج بازه ارتفاعی در هر زیرحوضه تعیین شد. نهایتاً ۳۹۶ واحد پاسخ 
هیدرولوژیکی در این مرحله ساخته شد. 
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پس از تعریف واحدهای پاسخ هیدرولوژیکی داده‌های هواشناسی 
به مدل معرفی گردید. در این مرحله پارامترهای هواشناسی مربوط به 
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در بازه زمانی ابتدای سال ۱۹۹۸ تا انتهای سال ۲۰۱۸ (مجموعاً ۲۱ 
سال) به مدل معرفی شد. همچنین آمار روزانه پارامترهای رطوبت 
نسبی, بارندگی» دمای حداقل و حداکثه ساعات آفتابی و سرعت باد 


نوری نژاد و همکاران. تحلیل حساسیت و عدم قطعیت مدل ٩۱۷۸‏ در شبیه‌سازی دبی جریان. رسوب و فسفر ... ۱۷۳ 


به‌عنوان پارامترهای هواشناسی داخل منطقه مورد مطالعه در بازه 
زمانی مشابه به مدل معرفی گردید. لازم به ذکر است تعدد داده‌های 
مفقود در اطلاعات هواشناسی به‌ویژه داده‌های روزانه دما قابل توجه 
بوده است. علاوه بر اطلاعات فوق. داده‌های دبی روزانه تخلیه شده از 
سد امیرکبیر و مشخصات فنی سد به مدل معرفی گردید. همچنین با 
استناد به اطلاعات اداره جهاد کشاورزی استان البرز اطلاعات 
مدیریتی اراضی باغی از جمله دوره‌های آبیاری و اطلاعات مربوط به 
کودهای فسفاته مصرفی در باغات منطقه نیز به مدل معرفی شد. پس 
از معرفی تمامی اطلاعات ورودی» جداول پایگاه داده مربوط به 
داده‌های ورودی مدل توسط مدل ساخته شد. در انتهاء اجرای نهایی 
مدل 5۱۷۸۲ صورت گرفت و نتایج شبیه‌سازی به‌صورت ماهانه و از 
ابتدای سال ۱۹۹۸ تا ۲۰۱۸ میلادی در پوشه خروجی پروژه ذخیره 
گردید. لازم به ذکر است که دو سال اول شبیه‌سازی به‌عنوان دوره 
آماده‌سازی" تعیین گردید. هدف از درنظرگرفتن این دوره دستیابی به 
ثبات مدل و حذف اثر داده‌های اندازه‌گیری مفقود و با خطادار می 


باشد. 


واسنجی, تحلیل عدم قطعیت و اعتبارسنجی 

ازآنجایی که با استفاده از مدل ۹۲۷۸۲ امکان تعریف و اعمال 
آبریز و نیز انجام برنامه‌ریزی درازمدت برای حوضه ایجاد می‌شود. 
نتایج باید قابل اطمینان بوده و دارای حداقل عدم قطعیت باشد 
(2007 ,.آه ۶1 ۲ا۸002:۲0). با توجه به متفاوت بودن بازه زمانی 
اطلاعات اندازه‌گیری شده متغیرهای دبی جریان» رسوب و فسفر کل 
در هر ایستگاه». حدود ۷۰ درصد مقادیر اندازه‌گیری شده برای دوره 
واسنجی و ۳۰ درصد آن برای دوره اعتبارسنجی مورد استفاده قرار 
گرفت که در جدول ۱ بازه‌های زمانی تعیین شده در هر ایستگاه 
مشخص شده است. به‌منظور واسنجی و تحلیل عدم قطعیت 
۶۷۷۵۸۲1۳ مورد استفاده قرار گرفت ( .۵1 6 تنامحعداخ۸ 
۵4, با استفاده از این الگوریتم عدم قطعیت خروجی توسط 
بازه‌های عدم قطعیت پیش‌بینی ۹۵ درصد» کمی‌سازی می‌گردد. در 
این روش با استفاده از فاکتور 18 که معرف پهنای بازه عدم قطمیست 
بوده و همچنین فاکتور ۳ که معرف نسبت مشارکت داده‌های 
اندازه‌گیری در بازه عدم قطعیت می‌باشد» کیفیت واسنجی مورد 
آرزیای قرار می‌گیرد.به این صورت که هرچه پهنای بازه عدم قعلمیت 
به عدد صفر و مقدار فاکتور ۴ به عدد ۱ نزدیک‌تر باشد. واسنجی از 
کیفیت مطلوب‌تری برخوردار است. برای مثال در مورد پارامتر دبی 


۱۳۹ 


فاکتور ۳ بزرگ‌تر از ۰/۷۰ يا ۰/۷۵ نشان‌دهنده کیفیت مطلوب 
واسنجی می‌باشد (2015 .4 6۶ تنامحعدحاحاظ) 

در نرم‌افزار ٩۷۷۸۲01/۳‏ برای شناسایی حساس‌ترین 
پارامترهای مدل به‌منظور استفاده در فرایند واسنجی دو روش تحلیل 
حساسیت چهانی و تخلیل حساسپت مطلق ‏ تعبیه شده استه در روش 
تحلیل حساسیت جهانی حساسیت پارامترها با استفاده از یک روش 
رگرسیون چندمتفیره محاسبه شده سپس از آزمون 1-7650 برای 
تعیین آاهمیت نسبی هر پارامتر استفاده می‌شود. در این روش 
حساسیت متغیر نسبت به تغییرات یک پارامتر درحالی‌ که باقی 
پارامترها نیز در حال تغییر هستند» سنحجیده شده و میزان حساسیت 
آن‌ها با ضرایب حساسیت 9]21-] و ۳-۷۵16 مشخص می‌گردد. هرچه 
قدر مطلق مقدار 121ه-] بزرگ‌تر و 0-۷72106 کوچک‌تر باشد پارامتر 
مورد نظر حساس‌تر است. درحال‌که در روش تحلیل حساسیت مطلق» 
حساسیت هر پارامتر نسبت به مقدار ثابت شده بقیه پارامترها سنجیده 
می‌شود (2007 ,.۵1 ۵ تتا0حواحض). تجربه ثابت کرده است که 
روش تحلیل حساسیت جهانی برای شناسایی حساس‌ترین پارامترها و 
تحلیل حساسیت مطلق برای بدست آوردن محدوده‌های مناسب 
پارامترها روند واسنجی مدل را تسهیل می‌کند ( ,.۵1 6۶ ۳۵۲2۳02721 
7 در این تحقیق با توجه به مرور مطالعات قبلی درباره 
واسنجیدبی جریان» رسوب و فسفر کسل در مدل 59۸ 
(2007 .6 6 تناموع9نظ :2010 .له ۶ 12160122060" : 
0 .۵ 6 حعفاعطامع؟ :2015 باه ۶ تنامجمودادض)» در 
مجموع ۴۰ پارامتر برای انجام تحلیل حساسیت جهانی در ابتدای 
فرایند واسنجی به‌منظور شناسایی حساس‌ترین پارامترها نسبت به سه 
متفیر مذکور و به کارگیری آن‌ها در واسنجی مدل تعیین شد. سپس 
برای تعیین بازه حساس هر پارامتر به‌منظور استفاده در فرایند 
(2007 ,.۵1 ۵1 تنامجعدداحه) مبنای کار قرار گرفت. 

در این مطالعه شاخص نش-ساتکلیف (۷5) مبنای اصلی ارزیابی 
کیفیت واسنجی مدل بوده و علاوه بر آن از ضریب تعیین (122 
مربعات نرمال شده (15۷۳) نیز استفاده شده است. از ۱۳۷۲5۳ 
در تحقیقات به‌عنوان نسبت انحراف میانگین مربعات خطا به انحراف 
استاندارد داده‌های اندازه‌گیری شده (56[۴) یاد شده است ( 1۷0۲1261 
7 , .01 61). 
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جدول ۱- بازه‌های زمانی تعیین شده برای دوره‌های واسنجی و اعتبارسنجی در ایستگاه‌های منطقه مورد مطالعه 
۵ 5000 ۲06 حرز عصمتاهاه 24 02410ناه۲۵ 20 ممتاعهانای [06مصه ۲0 )56 ۳6۵۲۵06 -1 12016 


دوره اعتبارسنجی دوره کالیبراسیون 
۱( ۳۵۱0 ما تاجن 
7 0 2012 1 0 2000 
6 0 2011 0 0 2000 


[111-2016 10 ۷۲2-7 


701-2016 0 66۵-8 


جدول ۲ طبقه‌بندی صحت نتایج مدل را بر اساس سه شاخص 
و ۳2 و 591 نشان داده که در ۳ کیفیت مدل‌سازی با چهار 
دسته نامطلوب» رضایت‌بخش, خوب و خیلی خوب تقسیم‌بندی 
می‌شود (2007 .]0 ۶۶ [2015(۷]01126 ,.۵1 6 ۷۲011251[ :1 


تحلیل حساسیت مدل 

کر انا یل سای خفانتی بهتطور شتاسایی بارافترهاین 
حساس مدل با توجه به بازه مقادیر مجاز تعریف شده در مدل 
۲ انجام شد. در شکل ۴ از میان ۴۰ پارامتری که برای رسوب 
و فسفر کل به‌صورت جداگانه ۱۷ پارامتری که قدرمطلق 21اه 
بزرگ‌تر و ۲-۷۵106 کوچک‌تری نسبت به سایر پارامترها داشتند و به 
عبارتی حساس‌تر بودند. نشان داده شده است که با درنظرگرفتن 


1۲60-2014 )0 10۵-6 


۷۲2۲-2014 0 10-6 


نام ایستگاه 
مصول ممتاداده 


متغیر اندازه‌گیری ده 
۵۸ ۱۵0وهع]۷[ 


دبی جریان و رسوب کیت 
عصمصصنع6 مه «ما۲ هعمج 
دبی جریان و رسوب سیرا-کرج 
۲صعصتله6 هه «ما۲ [ت-12٩‏ 

فسفر کل دیزین 

کناتمطام‌ومطظ 10121[ 1210[ 
فسفر کل شهرستانک 

قنامصامومطا م0[ کلقصهاده‌تطهط 


حساسیت بیشتر برخی از پارامترها به بیش از یک متغیر در مجموع 
۶ پارامتر با بیشترین حساسیت معرفی شد. بر اين اساس پارامترهای 
ضریب جریان پایه آب زیرزمینی (8۳ 4۸1۴11۸ ضریب زبری 
رواناب سطحی (0۷) و آهنگ افزایش بارش (۳۸۳۹) به 
ترتیب به‌عنوان حساس‌ترین پارامترها نسبت به دبی جریان» رسوب و 
فسفر کل شناخته شدند. بر اساس تحلیل حساسیت انجام شده 
پارامترهای آهنگ افزایش دما و بارش (۲۲۸۳5 و 5ضشا) و 
همچنین پارامترهای برفی مانند دمای بارش و فاکتور ذوب برف 
(5۷۳۳, ۹۳۲۷۲۳) جزء پارامترهای با حساسیت زیاد بوده‌اند که 
در مطالعات مشایهی که در مناطق کوهستانی انجام شده نیز جزء 
پارامترهای حساس به‌حساب آمده‌اند (2013 .01 6 حقحططفگ : 
0 , .۵1 6 هاا2ظ). 


جدول ۲- طبقه‌بندی کیفیت نتایج واسنجی بر اساس شاخص‌های آماری ۷5 ۳2 و 3616 (2007 ,۵ 61 نومنیم/(2015 ,.آه 1 نودنت0 :) 
تممتهم۷ :2007 ,۵ 6۶ تمدنم۷) وع‌نه‌صز ام‌ناعتاداه 1361۴ 0صد 1۳2 رک ده 0ععده مااباععد ممناهرطنلی ۵۶ ممنای‌تلتعععل بانامدم) -2 12016 


معیارهای ارزیابی عملکرد 
۱۹ 
رضایت‌بخش خوب 
ریت نان ۰0 


نامطلوب 


[۱0 ۲ 


0.70 > ۱5 > 0 0.50 > 5 > 0 ۱۲5 > 050 
0.70 > ۱5 > 0 0.45 > ۱5 > 0 ۱8 > 5 
0.50 > ۱۲5 > 5 0.35 > ۱۲ > 0 ۱۲5 > 5 
0.75 > ۳ > 5 0.60 > 8 > 5 1 > 0 
0.65 > 7 > 0 0.40 > ۳۳ > 5 1 > 0 
0.5 > 9۴ > 6 0.5 > 59۴ > 7 1۹1 < 7 


(2015 .۵1 1 
متیر معیار خطا 
خیلی خوب ۱۲۷6 دنت ۱ 
0 6۲۷ ۷ 
0 < ۱5 دبی جریان 
۳10۳ 
0 < ۱5 9 5 
احعصطنلع5 
65 < ۱۲5 #ر 
نموم 
5 < 82 دبی جریان 
۳10۳ 2 
0 < 132 ۳ 
مدمه ملع 5 
با ً با فسغة 
5 > 68 > 00 دبی جریان/ رسوب / فسفر 15۳ 


کمطام‌کمدام /احمحصتل۵۵/ ۳۱0۱۱ 
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نت 


ایا ۱ 


۷ 
۱۹ 


شکل ۴- نتایج تحلیل حساسیت پارامترهای مدل ٩۲۷۸۲‏ برای متغیرهای دبی جریان (ج)» رسوب (0) و فسفر کل 0( 
و) کتاا0دامو0 [062] 0ج (ظ) )ملع وره) ۲۱۵۲ ۵ ممامصصه ۳2 ۱۷۲۵۵۵ ۲۷۲۵۸۵۲ ۵۲ وتووافصق از لاتعصهو ۵۲ دااجمع1۳ -4 ۲۱۵0۲۵ 


این در حالی است که در مناطق غیر کوهستانی و پست. این 
پارامترها به‌ندرت به‌عنوان پارامترهای حساس محسوب شده‌اند 
(2013 .۵1 6 ۵1.۵ ۵ تتاو[ع۱۷]6۸ :2013 .۵1 6۶ ۱۱1۳2012 : 
9 در نتایج تحلیل حساسیت پارامترها برای متفیر فسفر کل نیز 
سه پارامتر ۳۸۳5 ,۲1۸۳5 و 8۳۳]1۷۳ از حساسیت بالایی 
برخوردار هستند. این موضوع اهمیت بالای دو عامل تغییرات ارتفاع و 
برف که از ویژگی‌های بارز مناطق کوهستانی است و نقش موثر آن‌ها 
در رژیم کمی و کیفی رودخانه در حوضه مورد مطالعه را نشان 
می‌دهد. 


مقادیر نهایی پارامترها در واسنجی مدل 

از مجموع ۲۶ پارامتر ذکر شده در شکل ۴ که حداقل نسبت به 
یکی از متغفیرها حساسیت قابل‌توجهی داشته‌انده ۲۳ پارامتر (جدول ۲) 
برای واسنجی مدل به کار گرفته شده و چهار پارامتر ضریب مانینگ 
بسرای کانال اصلی (۸2 013 فساکتور فرسایش کانسال 


(013_۳0۳(/0) شاخص دسترس‌پذیری فسفر (۳۹۴) و غلظت 
رسوب در جری‌ان جانبی و آب زیرزمینی (1_8ه) با 
درنظرگرفتن حساسیت کمتر آن‌ها در مقایسه با ۲۲ پارامتر دیگر در 
واسنجی مدل مورد استفاده قرار نگرفتند. مقادیر نهایی پارامترها در 
واسنجی پس از انجام ۱۷ تکرار در ٩۷۷۵۸۲0/۳‏ بدست آمد. 
مطابق جدول ۴ مقادیر پارامترهای مربوط به برف و آهنگ افزایش 
دما و بارش (شماره ۱ تا ۵) پس از تکرار پنجم در فرایند واسنجی و 
دستیابی به مقادیر نسبتاً مناسب شاخص‌های آماری (جدول ۲) بر 
اساس مقادیر نمونه گیری شده از بهترین دور شییه‌سازی مدل ابت 
گردید. این عمل باعث کاهش قابل توجه ضخامت بازه عدم قطعیت 
۵ درصد (فاکتور 1۶) برای هر سه متغیر و نزدیک شدن به مقدار ۱/۵ 
و کمتر از آن شد که برای مدل‌سازی مطلوب توصیه شده است 
(2015 ,.۵1 2 نعد1ت۷]0). سپس به‌منظور دستیابی به نتایج بهتر با 
ییورارههای #دررق قنه برای سای پا اسرد تگرارهای پمایفز 
نهایت حدود نهایی آن‌ها در آخرین تکرار بدست آمد. 
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جدول ۲- مقادیر واسنجی شده پارامترها 
۲۷۵۵66 021۵۵۵۲ 021101۵60 -3 م1201 


روش تغییر 


و ۷( 0 ِ 
9 
[0.20,0.35] ۱ اع. 02 12 
[0.09,0.12] 2-۵۵ بط و ط اما 13 
[0.79-,0.85-] 7-۵۵ 19_0۱ 14 
[0.10,0.18] 7-۵۵ ۱( 15 
[3,5] 6+( 066 16 
[0.01-,0.12-] 7-۵۵ ۹ 0 0 ((0] 17 
[0.53,0.57] ۱۹ ۲۵ 0۳ 18 
[0.84,0.94] ۱۹ 00 0۶ 19 
[0.007,0.008] 6+( مد 5۳0۵ 20 
[1.33,:1.39] 2-6( موه 9۳۳۲ 21 
[23,0.61 .0] ۱۹ 10 22 
ارزیابی واسنجی دبی جریان 


بر اساس شکل‌های ۸۵ ۶ و ۷ در سه ایستگاه هیدرومتری مورد 
مطالعه مقدار 120107- کمتر از ۱/۵ و مقدار 2010۲/-۳ بیشتر از ۰/۷ 
بدست آمده است که بیانگر کیفیت مطلوب واسنجی دبی جریان است 
(2015 ,.۵1 ۶ ۲اهعوح00ض۸). با توجه به جدول ۴ در دوره واسنجی و 
اعتبارسنجی مدل و در اغلب موارد 5 و "15 به ترتیب بزرگ‌تر از 
۵ ۰/۶ بوده و با ۷۲9۳( بین ۲/۲۴ تا ۸/۸۲ مترمکعب در ثانیه 
مقدار 518 کمتر از ۰/۷ به دست آمده که بر اساس طبقه‌بندی کیفی 
مدل‌سازی (جدول ۳) مدل عملکرد رضایت‌بخش و خوب داشته است 
(2015 ,اه ۶ ۷022). در ایستگاه بیلقان و به‌ویژه در دوره 
واسنجی مقادیر این شاخص‌ها نشان از پیش‌بینی ضعیف‌تر دبی جریان 
به نسبت سایر ایستگاه‌های مشاهداتی دارد. با توجه به مقایسه مقدار 
فاکتور *1 در شکل‌های ۸۵ ۶ و ۷ خروجی‌های مدل در ایستگاه بیلقان 
دارای عدم قطعیت بیشتری نسبت به سایر ایستگاه‌های سنجش دبی 
هستند. بر اساس جدول ۵ که متوسط متغیرهای مشاهداتی و 
شبیه‌سازی شده در فصول مختلف سال را نشان می‌دهد. به جز فصل 
تابستان» در باقی فصول مقدار دبی کمتر از مقدار واقمی پیش‌بینی 
شده است. همچنین بیشترین تفاضل بین مقادیر مشاهداتی 9 
شبیه‌سازی شده مربوط یه فصل بهار می‌باشد که با توجه به شکل‌های 
۵ ۶و ۷ در نقاط اوج در ماه مه میلادی (اواسط بهار) این اختلاف 
قابل ملاحظه است؛ بنابراین عمده ضعف مدل در هر سه ایستگاه عدم 
پیش‌بینی درست نقاط اوج جریان بوده است. این مسئله را می‌توان به 


روش تغییر 
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2۵۶ ۲2 معصمطل 1( .0( 
4.44 تاره ۱۵۵ مو .5۳۲/۲۳ 1 
094 تاره ۱۵۵ صو.1۷]۳۴ 1[ 2 
060 تاه ۱۵۵ ۷۳۷/۸ 3 
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[173,329] ۱۹۵ 7 0 7 ( (6 7 
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[0.11,0.32] ۱۹ ۷ ۲۲۸ ,]۸ 9 
۲197:2701 ۱۹۵ 0۳16 10 
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عم پیش پیتی صاخ تربار ورودت ب رردحاخهردر آکر دوب ره 
تعمیم داد که تقریبا از اواسط بهمن (فوریه میلادی) روند افزایشی آن 
شروع شده و اواسط اردیبهشت (مه میلادی) به اوج خود می‌رسد. چند 
عامل می‌تواند در عدم تشخیص صحیح این نقاط موثر باشد. اول 
اینکه تعدد قابل توجه داده مفقود در ایستگاه‌های هواشناسی می‌تواند 
بر عدم پیش‌بینی صحیح مدل موثر باشد ( ,.۵1 6 طدنصدادم؟۳ 
05 که در اين منطقه بخصوص کمبود داده‌های دمای روزانه در 
مقایسه با سایر پارامترهای هواشناسی قابل توجه بوده است. دوم اینکه 
مقادیر نهایی مربوط به پارامترهای دمای بارش و ذوب برف و آهنگ 
افزایش دما و بارش به کل زیرحوضه‌ها تعمیم داده شده و احتمالا 
لازم است تا مقادیر آن‌ها با توجه به عدم قطعیت زیاد این پارامترها در 
تواحی کوهستانی (2010 .۵ تحتف!1ن ,همعط 0صه متلف؟ : 
6 ,00۵۱ 0طه حعص :2011) در سطح زیرحوضه‌ها تنظیم 
شود. اما بهترین پیش‌بینی‌های مدل با توجه به جدول ۵ در فصل پاییز 
بوده که اختلاف بین مقادیر شبیه‌سازی شده و مشاهداتی کمتر از 
مقدار متناظر برای سایر فصول بوده است. این نشان می‌دهد که 
پارامترهای کنترل کننده جریان پایه از جمله ضریب جربان پایه آب 
زیرزمینی (11 ۸1۳۲۱۸ زمان‌ت-أخیر آب زیرزمینی 
(0۳۲۸۷ 0۷۷ عمق آستانه آب در آبخضوان کم‌عمسق 
(0۷۷۵((۷) و عمق آستانه آب برای حرکت از سفره آب کم‌عمق به 
آب‌های زیرزمینی (181۷۸۳۷۲) که جزو پارامترهای حساس مدل 


هستند به‌خوبی تخمین زده شده‌اند. 


نوری نژاد و همکاران. تحلیل حساسیت و عدم قطعیت مدل ٩۷۷۸1‏ در شبیه‌سازی دبی جریان. رسوب و فسفر ... 
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۱۷۸ آ 


نشریه آب و خاک. جلد ۰۳۶ شماره ۰۲ خرداه - 


تیر 


۱۴۳۰ 


نوری نژاد و همکاران. تحلیل حساسیت و عدم قطعیت مدل ٩۱۷۸‏ در شبیه‌سازی دبی جریان. رسوب و فسفر ... ۱۷٩‏ 


رضایت‌بخش و خوب داشته است. با توجه به جدول ۵ مدل در اکثر 
شهرستانک» مقدار رسوب را کمتر از مقدار واقعی پیش‌بینی کرده 
است. مانند دبی» بیشترین اختلاف بین مقادیر مشاهداتی و شبیه‌سازی 
اوج رسوب واضح‌تر به نمایش درآمده است. بهترین پیش‌بینی‌ها نیز 
مربوط به فصل پاییز است که داری کمترین مقدار رسوب در میان 
تمام فصول بوده و تفاضل بین مقادیر مشاهداتی و شبیه‌سازی شده در 
کمترین حالت قرار دارد. مقایسه شاخص‌های آماری به دست آمده 
به ویژه ولا برای متغیرها نشان می‌دهد که شبیه سازی رسوب سبت 
به دبی جریان با دقت کمتری انجام شده که با مطالعات مشابه در 
زمینه واسنجی دبی جریان و رسوب همخوانی دارد ( 200 ۳۲۱۶ 
1 ,167۷ ص2 2019۷ رصوت۷ 20 معد1210«7 :). دلیل این 
موضوع می‌تواند این باشد که با توجه به تأثیرپذیری رسوب از رژیم 
جریان (2017 ,.۵1 61 1<۸۲0510) نه‌تنها عوامل موّثر در کاهش 


13-۳ 20101 2 4 

01 < ۳2610۲- 
اعتبارسنجی واسنجی 
۱ 4 «متاحتطانلدت 


سم پسي هسي پسي پسي سم هپس پستي سیم نس پم 


سری زمانی 
2۵ 11۳6 متاهداتی لعبعو00 


دقت پیش‌بینی جریان می‌تواند در تعیین میزان بار رسوب منتقل شده 
موّثر باشد» بلکه وابستگی مدل ٩۷۷۸۲‏ به بسیاری از مدل‌های 
تجربی و نیمه‌تجربی مانند ۹05-0 و تاون]۷ را می‌توان عامل 
موّثری دانست که باعث می‌شود تا ٩۷۷۸۲‏ بار رسوب رابادقت 
کمتری ردیابی کند (1970 ,۹0116 20 طوعآ). علاوه بر این از 
آنجایی که مقادیر تخصیص داده شده به پارامترهای حساس مرتبط 
با روندیابی رسوب و فرسایش مانند پارامتر خطی محاسبه‌کننده 
حداکثر رسوب مجدداً وارد شده به کانال ([۳00) فاکتور توان در 
روندیابی رسوب در آبراهه (5۳17۳) فاکتور فرسایش‌پذیری خاک در 
معادله ]79 (1_آ9لا)» فاکتور حفاظت خاک در معادله ]نا 
(۳_طا5نا)» فاکتور فرسایش کانال (0۳1-00۷1) و همچنین سایر 
پارامترهای حساس از جمله ضریب زبری رواناب سطحی (0۷-3) 
به کل حوضه تعمیم داده شده است به نظر می‌رسد با توجه به تنوع 
زیاد شیب و ارتفاعات در منطقه» ارزیابی و واسنجی این پارامترها در 
سطح زیرحوضه بتواند عملکرد مدل را بهبود ببخشد. 
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زه عدم قطعیت ٩۵‏ درصد 9500 


مشاهداتی 00:6۳:20 


شبیسازی شده 64)قاناص۹1 


شکل -٩‏ نتایج واسنجی و اعتبارسنجی رسوب در ایستگاه سیرا - کرج 
ماه [51۳۵-1622 26 و۲6۶6 60 ۷۵۱۱02 0صره مت موانلی )حمصنلع6 -9 ۲12۱۲۵ 


۰ نشریه آب و خاک. جلد ۰۳۶ شماره ۲. خرداد - تیر ۱۴۰۱ 


ارزیابی واسنجی فسفر کل در فصل زمستان بیشتر از آن برآورد گردیده 

فاکتور 1 در هر دو ایستگاه دیزین و شهرستانک کمتر از مقدار 
متناظر آن در پیش‌بینی دبی جریان و رسوب است که نشان‌دهنده 
عدم قطعیت کم پارامترهای استفاده شده در واسنجی فسفر کل است 
وه دلیل فقدان کمتر فاکتور 1 مسبت متا کت ففادیر مشباهداتن داشته و بیشترین اختلاف مقادیر مشاهداتی 
فسفر کل در بازه عدم قطعیت ۹۵ درصد (فاکتور ) نیز عددی 


بین مقادیر مشاهداتی ۳ شبیه‌سازی شده برای ایستگاه دیزین در فصل 
زمستان و برای ایستگاه شهرستانک در فصل پاییز رقم خورده است 
که ضخامت بازه عدم قطعیت ٩۵‏ درصد در کمترین حالت خود قرار 


وان شبیه‌سازی شده نیز 


مربوط به فصل بهار بوده که در نقاط اوج (مه میلادی) این اختلاف به 


کوچک‌تر به دست آمده است. یکی از دلایل کمتر بودن فاکتور ۲ در بیشترین حد خود می‌رسد (شکل‌های ۱۰ و ۱۱). فسفر کل در این 


پیش‌بینی فسفر کل» پهنای بسیار کم بازه عدم قطعیت (فاکتور کل در 


منطقه با توجه به حاصل شدن ۱5 بالای ۰/۹ بسیار نزدیک و یا بهتر 


ماه‌هایی است که مقادیر فسفر کل کمتر از ۲۰ تن می‌باشد (شکل‌های از برخی مطالعاتی بوده که در حوضه‌های کوهستانی یا غیر کوهستانی 


۱ ۱۱). با توجه به نتایج ارائه شده در جدول ۴ که شاخص‌های 
5 و ۳ بیش از ۰/۹ بوده و با ]۷ بین ۱/۲۸ تا ۲/۷۴ تن در 


مدل‌سازی کرده‌اند (2013 ,.۵1 ۵ قانهتذل :2016 ,.۵1 61 بان : 
6 راتدطاها۳۲۵ 200 6۲ع261). این نشان می‌دهد که ٩۷۷۸۲‏ در 


ماه حقدار 6۳ در نو استگاه بن ۰/۱۷ ۵ ۳۵/عاصهان شنده. اقا شبیه‌سازی فسفر کل در مناطق کوهستانی مشابه هم از نظر 
۳ از 97 یبن 


در پیش‌بینی فسفر کل عملکرد بسیار خوبی داشته و بهتر از دبی 
جربان و رسوب تخمین زده شده است. با توجه به جدول ۵ فسفر کل در رزوی مت نیز 


به طور متوسط در اغلب فصول کمتر از مقدار واقعی پیش‌بینی شده اما 


زمان‌بندی تغییرات بار فسفر کل و هم از نظر مقدار تغییرات آن 
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شبیه‌سازی شده 4عاهاراحط۹ 


شکل ۱۱- نتایج واسنجی و اعتبارسنجی فسفر کل در ایستگاه شهرستانک 


صمناهای عاقصدای طاحطاه اه عاایاعع۲ «مت)هنا۷۵ 0صه ممتاهرطانای کامطام‌عمدام 10621 -11 م۲1 


نوری نژاد و همکاران. تحلیل حساسیت و عدم قطعیت مدل ٩۱۷۸‏ در شبیه‌سازی دبی جریان. رسوب و فسفر  ...‏ ۱۸۱۲ 


نتیجه گیری 

بر اساس نتایج تحلیل حساسیت. پارامترهای ضریب جریان پایه 
(0۷-۱۷) و آهنگ افزایش بارش (۳],۸۳5) به ترتیب حساس‌ترین 
پارامترها نسبت به دبی جریان, رسوب و فسفر کل بودند. همچنین 
در تحلیل حساسیت مدل برای متفیر فسفر کل» سه پارامتر 
۳۳۸5 ,۲1۸۳5 و ٩۳۲/۴‏ از حساسیت بالایی برخوردار بودند. 
این موضوع اهمیت بالای دو عامل تغییرات ارتفاع و برف و نقش 
موثر آن‌ها در رژیم کمی و کیفی رودخانه در حوضه کوهستانی مورد 
مطالعه را نان می‌دهد. با توجه به ضرایب ۵15 *12 3۷6۳ و 
5۴ بدست آمده در دوره واسنجی و اعتبارسنجی هر سه متغیر 
رواناب» رسوب و فسفر کل» مدل در شبیه‌سازی رواناب و فسفر کل 
عملکرد بهتری داشته است. به‌ویژه در مورد فسفر کل و در دوره 
واسنجی مدل ضرایب ولا 9 ۳ دارای مقادیر بالای ۰/۹ و 4 
کمتر از ۰/۵ بوده که نشان‌دهنده پیش‌بینی بسیار خوب مدل است. 
با توجه به مشارکت بالای ۷۰ درصد مقادیر اندازه‌گیری شده رواناب 
و رسوب و حدود ۶۰ درصد مقادیر اندازه‌ گیری شده فسفر کل بازه 
عدم قطعیت منطقی و مناسبی برای خروجی‌های مدل به دست آمده 
که ناشی از انتخاب محدوده مقادیر مناسب برای پارامترها در فرایند 
واسنحی است. در هر سه متغیر بیشترین عدم قطعیت مربوط به 
فصل بهار ب‌ویژه در نقاط اوج جریان (مه میلادی) بوده که شامل 
بیشترین مقادیر مشاهداتی در هر سه متغیر بوده و دقت پیش‌بینی 
مدل نیز کمتر می‌باشد. بهترین پیش‌بینی‌های مدل نیز عمدتا در 
فصولی است که دبی متوسط جریان به دبی پایه نزدیک بوده و 
کمترین بار ماهانه رسوب و فسفر در رودخانه منتقل می‌شود. با 
درنظرگرفین تفاضیل مقادیر متفسط مخاهداتی و شببه‌سازی شنده 
(جدول ۵) برای دبی جریان و رسوب در فصل پاییز و برای فسفر کل 
در فصل زمستان بهترین پیش‌بینی‌ها انجام شده که دارای کمترین 


منابع 


عدم قطمیت به نسبت فصول دیگر سال هستند. مشکل پیش‌پینی 
نقاط اوج دبی جربان و رسوب در مطالعات مشابه مربوط به 
شبیه‌سازی حوضه‌های آبریز کوهستانی در ایران و خارج از ایران نیز 
تجربه شده است (2008 .اه 6 صفتصصعاومک؟ ص۲۷2 4صح صصیراظ : 
90 ,۷29 هه 12107516 :2011 ,#«عنآ]) که می‌توان آن ر 
به عوامل مختلفی از جمله کمبود داده‌های دمای روزانه» وابستگی 
۲ به مدل‌های تجربی روندیابی رسوب و فرسایش, و عدم 
بررسی تغییرات برخی پارامترهای دارای عدم قطعیت به تفکیک 
زیرحوضه‌ها در واسنجی مدل مرتبط دانست. 

با توجه به اينکه اکثر شاخص‌های ارزیابی کیفیت» حاکی از 
عملکرد مطلوب شبیه‌سازی دبی جریان» رسوب و فسفر کل در ایین 
حوضه بوده‌اند» مدل واسنجی شده در این تحقیق می‌تواند مبنای 
پیش‌بینی‌های آتی در خصوص کمیت‌های مورد اشاره قرار گیرد و بر 
اساس آن می‌توان سناریوهای مدیریتی برای کنترل کیفیت آب 
تعریف کرد. از جمله محدودیت‌ها در تحقیق حاضر می‌توان به وجود 
تعداد زیاد داده‌های مفقود در ایستگاه‌های هواشناسی و بخصوص 
برای متغیر دما اشاره کرد که به نظر می‌رسد با توجه به تغییرات زیاد 
دما و بارش در ارتفاعات به دلیل ماهیت کوهستانی منطقه» هرچه 
اطلاعات کامل‌تری از اندازه‌گیری‌های روزانه ایستگاه‌های هواشناسی 
به مدل ارائه شود. مدل پیش‌بینی بهتری از روند تغییرات کمیت‌ها 


سپاسگزاری 
نویسندگان این تحقیق از سازمان حفاظت محیطزیست اداره 


می‌آورند. 
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